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ABSTRACT 

We investigate possible environmental and morphological trends in the z ~ bar fraction using two 
carefully selected samples representative of a low-density environment (the isolated galaxies from the 
AMIGA sample) and of a dense environment (galaxies in the Virgo cluster) . Galaxies span a stellar 
mass range from 10^ to lO^^Mo and are visually classified using both high-resolution NIR (H-band) 
imaging and optical rgb images. We find that the bar fraction in disk galaxies is independent of 
environment suggesting that bar formation may occur prior to the formation of galaxy clusters. The 
bar fraction in early type spirals {Sa — Sb) is ~50%, which is twice as high as the late type spirals 
{Sbc ~ Sm). The higher bar fraction in early type spirals may be due to the fact that a significant 
fraction of their bulges are pseudo-bulges which form via the buckling instability of a bar. i.e. a 
large part of the Hubble sequence is due to secular processes which move disc galaxies from late to 
early types. There is a hint of a higher bar fraction with higher stellar masses which may be due 
to the susceptibility to bar instabilities as the baryon fractions increase in halos of larger masses. 
Overall, the SQ population has a lower bar fraction than the Sa — Sb galaxies and their barred 
fraction drops significantly with decreasing stellar mass. This supports the notion that SOs form via 
the transformation of disk galaxies that enter the cluster environment. The gravitational harassment 
thickens the stellar disks, wiping out spiral patterns and eventually erasing the bar - a process that is 
more effective at lower galaxy masses. 
Subject headings: galaxies:barred galaxies 



1. INTRODUCTION 

Understanding bar-formation and the secular effects 
that bars have on the stellar component is becoming cen- 
tral to our understanding of galaxy formation and evo- 
lution. Once a bar forms it can change the scale length 
of the stellar component via scattering/mixing of stars 
in the radial direction a nd they can create an extended 
stellar distribution (Roskar ct al. 2008). Bars can move a 
galaxy between morphologic al classes through the secular 
formation of pse udo-bulges ( Athanassoula] 2002 [Debat- 
tista et al. 2004 ) and they can drive gas to the central 
black hole fuelling AGN (Shlosman et al. 1989). 

Numerical simulations have shown that ba rs natu- 
rally arise from the secular evol ution of discs (iT^oomre 
19641 lOstriker fc Pe ebles 1973; Fall fc Efstathiou||iy8U| ). 
Bars can also be tr iggered b y dynamical interacti ons in 
the field (|Gerin et al. 199 0j jBarnes fc He rnquist ]|1992 
Romano-Diaz et al.||2008| |Dubinski et al.||2009| ). It has 
been shown that close galaxy companions are associated 
with bar formation, but p rimarily for early Hubble-types 
( Elmegreen et al.]|l990| ). In galaxy clusters, gravita- 



stellar bar which is subsequently heated, cau sing the rem 



tional encounters (harassment) can drive a morphological 
transformation from late type disks to dwarf ellipticals 
(dSphs). In this scenario, encounters create a "naked" 
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nants to become more spherical with time (Moore et al. 
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Is environment the key factor in determining why two 
similar galaxies may or may not have a bar or is the ex- 
istence of a bar related to the initial conditions of galaxy 
formation? The stability (or instability) of disks to bar 
formation may also depend on the baryon fraction, and 
in particular the mass of stars and gas in the disk. This 
varies across t he Hubble sequence and depends strongly 
on halo mass (v an deri Bosch et al.|2000 Courteau et al. 



2003| |McGaugEfp005p . I'his all suggests that in addi 
tion to local density, morphology and halo mass are also 
important principal parameters to investigate. Another 
key observational result is the fact that the bar fraction 



is not changing significantly with red shift (se elE lmegreen 
et al |2004l]Mariiiova fc Jogee|2007l but alsolSlieth et aT 



2008 for a ditt'erent result), however, most disk galaxies 



are not within dense environments so it would be diffi- 
cult to disentangle the effects of environment, especially 
at higher redshifts. 

Several studies have shown that there is no evidence 
for a dependence of ba r frequency on galaxy environment 



van den Bergh 2002), the same is true even if galaxies 
of diffe rent morphol ogical typ e are cons i dered indepen- 



dently ( |Aguerri et al. 2009). |Li et al.| ( |2009| ) came to 
the same conclusion analyzing the clustering properties 
of barred and unbarred galaxies of similar stellar mass 
and finding it indistinguishable over all the scales probed 
(from ~20 kpc to 30 Mpc). More recently, the Coma 
cluster was studied by Mendez-Abreu et al. (20101: they 
find that the bar fraction does not vary significantly even 
when going from the center to the cluster outskirts. How- 
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Fig. 1. — H-band luminosity distribution of the Field and ViRGO 
moderately inclined (axis ratio > 0.4), morphologically selected 
discs (with Hubble type ranging from SO to Sm). Shaded his- 
tograms represent the barred discs. 



ever, the Coma cluster is such an extreme environment 
that most of its apparent spiral galaxy population may 
be field galaxies in projection. 

In the light of these observational results and motiva- 
tion from numerical simulation studies, we aim at mea- 
suring the bar fraction (as number of barred discs over 
the total number of discs) as a function of environment 
and disc morphology, at z ^ in two carefully selected 
samples representative of a low-density environment (the 
isolated galaxies from the AMIGA sample) and of a mod- 
erately dense environment (galaxies in the Virgo cluster). 

To achieve this goal it is important to use homoge- 
neous classifications since, as we have shown in Giordano 
et al., (2010) (paperl hereafter), the bar fraction is very 
stable against sample selection but that some (possibly 
spurious) differences can arise if the comparison is based 
on samples classified using different methods (for exam- 
ple visual classification versus automated profile fitting) . 
In particular, the way the disc population is identified 
and isolated plays a crucial role, since, if no detailed 
morphological information is available, discs can easily 
be miscounted (for example applying only color and/or 
magnitude cuts). 
In order to address this, we use data from the UKIDSS 



Large Area Survey (Lawrence et al. 2007 1 and from SDSS 
DR7 (Abazajian et al. 2009J , with the great advantage 
of combining optical rgb images with near-infrared (H- 
band) imaging with excellent resolution for local universe 
studies, that allow us to visually inspect the images to 
provide detailed morphological classifications. 

The outline of the paper is the following: in section S|2] 
we present the data that we are using, their classification 
and the selection of the samples based on local density 
estimation. The results about the bar fraction in the 
different cases are presented in section fjSjand discussed 
in section ^ 

2. DATA 

2.1. Virgo sample 

In Paperl we presented a thorough study of the barred 
galaxies in the Virgo Gluster from which we adopt all the 
classified galaxies with a measured H-band magnitude 
from 2MASS. The ViRGO disk sample is composed of 
moderately inclined (axis ratio larger than 0.4) members 
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Fig. 2. — Bar fraction as a function of local density measured via 
the p5 proxy of the NAIO (grey diamonds), Field (blue square) 
and Virgo (red circle) samples, with Field and and ViRGO samples 
bracketing the extremes in ps of the NAIO density distribution. 
Averaged over all morphological types, the bar fraction is constant 
against local density. 



with UKIDSS near-IR imaging of Hubble type between 
SO and Sm, spanning a H-band magnitude (stellar mass) 
range of -17 to -25 mag (10® to lO^^ y^^y ^j^e follow- 
ing analysis, we use the H-band magnitudes from Paper I 
to compute stellar masses assuming a flat {B — H) color 
with a Th,* = 1. The local gal axy density for mem bers 



is determined via the proxy (Baldry et al 
ing the p ositions and magnitu des from the Virgi 
Catalog ( [Binggeh et 'ar]|1984[ ) 



2006] ), us- 
;o Cluster 



2.2. Field sample 

To provide a robust comparison to the ViRGO sample, 
we select a true Field sample using the AMIGA (Anal- 
ysis of the interstellar Medium of Isol ated GAlaxies) 

The AMIGA 



2005) 



project (Verdes-Montenegro et al. 
catalo gue is based on the KIG catalog Karachentseva| 
(|1973| of isolated galaxies (z < 0.1). The KIG cat- 



alog is composed of 1050 galaxies with apparent blue 
magnitudes brighter than 15.7 mag; these isolated galax- 
ies are selected to have no neighbor of comparable size 
within twenty galactic diameters. The KIG catalog has 
been used by multiple studies to investigate the effects 



of under-dense environment on galaxy p roperties ( Adams 
et al.||1980[ [Haynes fc Giovanelli| ll980) and the AMIGA 



project quantified the isolation of KIG galaxies i dentify- 
ing th e ir samp le of 791 genuinely isolated galaxies | Verley| 
et al.| ( |2007a ). The AMIGA project has also compiled 
multiwavelength coverage of this statistically significant 
sample of the most isolated galaxies in the local universe 
and the dataset includes optical photometry and mor- 
phol ogies for a r edshift-complete subsample of 956 galax- 
ies (S ulentic et al.„2006) , and these data are pubhcly re- 
leased under a VO interface at http://amiga.iaa.es/. 

By cross-matching the AMIGA catalogue with the 
2MASS and SDSS databases, we identify 563 galaxies 
with both H-band and rgb images. As in our ViRGO 
sample, we select only moderately inclined (axis ratio 
> 0.4) disk galaxies with Hubble type of 50 to Sm. The 
stellar masses of the resulting 390 Field disk galaxies 
are determined using the total H-band magnitudes from 
2MASS and assuming the same flat {B — H) color and 
Th,* = 1 as in the ViRGO sample. The local density for 
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each gal axy is determi ned by their 5th nearest neighbors 
{P5, see [Baldry et al.| [2006') as defined by the AMIGA 
catalogue of neighbors, constructed down to a magnitude 
m ^17.5, lying with in 0.5 Mpc around the KIG galaxies 
( | Verley et aL]|2007b[ ). 



2.3. Identifying Barred Disks 

The UKIDSS Large Area Survey is an ongoing sur- 
vey to image 4000 deg^ at high Galactic latitudes in the 
YJHK filters to a depth in H of 18.8 mag; it has a spa- 
tial resolution of 0.4"/pixel (like SDSS data) and an av- 
erage seeing of 0.8". Our galaxies span a redshift out to 
z ~0.03 where the UKIDSS imaging has a physical res- 
olution of ^400 pc/pixel. As outlined in Paperl, we use 
the H-band imaging to visually classify the disk galaxies 
into one of three categories: "barred" , "non-barred" , or 
"uncertain" . All the galaxies in the ViRGO sample have 
H-band imaging from UKIDSS that is also available for 
172 galaxies in the Field sample, for the rest we must 
rely on SDSS z' and rgb imaging; however, for all the 
galaxies with both H-band and z' imaging, we find that 
our bar classifications are essentially identical. 

Both Virgo and Field catalogues, comprehensive of 
rgb and H-band thumbnails are available online^. 



3. RESULTS AND DISCUSSION 
3.1. Barred Fraction vs. Environment 

In Figure [l] we show the H-band luminosity distribu- 
tion for moderately inclined disk galaxies in our Field 
and Virgo samples. The barred fraction averaged over 
the disk population in the Field is --34% (133/390) and 
in Virgo is ~28% (90/332). We find that even with bar 
classifications based on high resolution near-IR imaging, 
the barred fraction does not vary with environment when 
considering the disk population as a whole. Although 
our results are consistent with earlier wor k based on op- 
tical imagin g (see for example van den Bergh 2002 Li 
et al. 2009), the constancy of the barred traction con- 
tlicts with expectations from current galaxy formation 
models, since strong interacti ons can trigger bar insta- 
bilities ( [Berentzen et al.|[2004l ). 

To further test the robustne ss o f our result, we in- 
corporate the Nair & Abraham (2010) optically-selected 
sample of about 14,000 galaxies from the SDSS DR4; 
this sample includes nearly all spectroscopically-targeted 
galaxies in the redshift range 0.01 < z < 0.1 to an appar- 
ent extinction-corrected limit of g < 16 mag. In addition 
to visually classified Hubble types, the NAIO catalogue 
contains the existence of bars , stella r masses computed 
according to KaufFmann et al. (2003), and an estimate of 
the local density comp uted~iising the 5th n earest neigh- 
bors (ps) according to Baldry et al. ( |2006 ). To ensure 
that we can compare the (optical) JNAIO sample directly 
to our (near-IR) results, we apply the same disk selection 
criteria as in our Field and ViRGO samples and consider 
only galaxies at z < 0.03. 

Figure [2] shows the barred fraction in our Field (blue) 
and Virgo (red) populations as well as the NAIO sample 
(grey) as a function of local galaxy density (^5); note 



^ http://www.itp.uzh.ch/ giordano/mwpbg-v2.htmI 



how effectively our Field and Virgo samples bracket 
the extremes in p^. We find the barred disk fraction is 
surprisingly resilient (~ 30 — 35%): the barred fraction 
does not vary as the environment changes from isolated 
field galaxies to cluster cores. 

3.2. Barred Fraction vs. Disk Morphology 

To investigate how the barred fraction varies with 
disk morphology, we divide both our Field and ViRGO 
samples into three classes: 1) lenticulars (featureless 
discs, corresponding to Hubble- types S'O — S'Oa); 2) early- 
type spirals (bulge-dominated discs, corresponding to 
Hubble- types Sa — Sb); and 3) late- type spirals (disc- 
dominated or bulge-less discs, corresponding to Hubble- 
types Sbc— Sm). In Figurejsjwe show the barred fraction 
as a function of stellar mass for the three disk classes 
where the blue squares represent the Field and the red 
circles the ViRGO members. Each panel also includes the 
average barred fraction for the three disk classes in the 
Field (dotted fine) and in ViRGO (dashed line). The 
differences in the relative number of galaxies in each disk 
class is due to the morphology-density relation, i.e. the 
fraction of lenticulars in the Field (< 10%) is lower than 
in Virgo (- 30%; Paperl). 

We find that the barred fraction for early-type spirals is 
systematically higher than in late-type spirals regardless 
of environment: 45 — 50% for Sa — Sb vs. < 25% for 
Sbc — Sm. The barred fraction in lenticulars (50 — S'Oa) 
is also lower in both environments. 

4. SUMMARY 

We present the first comprehensive study of barred 
disks as a function of environment that uses NIR and 
rgb imaging to resolve bars; the advantage of using near- 
infrared imaging from UKIDSS is that bar classifications 
are less affected by dust and bright star-forming regions. 
We expand on our study of bars in ViRGO by building a 
Field sample using the KIG catalog of isolated galaxies. 
Our Field and Virgo disk populations are at z < 0.03, 
span a range in stellar mass from ^lO^to^lO'^^ Mq and 
Hubble type (SO — Sm), encompass a wide range in local 
densities and are analyzed in exactly the same manner. 

We find that the barred disk fraction is surprisingly 
constant at 30 — 35% in both the Field and Virgo sam- 
ples, i.e. the barred fraction for the disk population as 
a whole does not depend on environment. We test the 
robustness of our result by analyzing the NAIO optically- 
selected sample of nearby galaxies in the same manner, 
and we again find a constant barred fraction across the 
full range of local galaxy density. 

This implies that disks become barred prior to the 
late time assembly of galaxy clusters, which is consis- 
tent with observational evidence that the bar fraction 
does not evolve strongly with redshift. 

The barred fraction is highest for early-type spirals 
{Sa — Sb) regardless of environment: these galaxies are 
nearly twice as likely to be barred as late-type spirals 
{Sbc — Sm). If a late type spiral forms a bar, then it 
may also form a pseudo-bulge via a buckling instability 
and its morphological class will change. Indeed, the con- 
sensus is forming that our own Galax y has evolved acros s 
the Hubble sequence in this fashion ( Agertz et al.|20l"0 ). 
If this is a common ph enomenon , as numerical simula- 
tions indicate (Debattista et al.|2006^ , then we naturally 
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Fig. 3. — Bar fraction vs. stellar mass distribution of the Field and ViRGO morphologically selected disks. The red circles represent 
the Virgo sample, the blue squares the Field sample. In each panel, the average value over the morphological type is represented by the 
dashed/dotted line, for the ViRGO (red) and the Field (blue) respectively. The Sa — Sb disks have the highest bar fraction, regardless of 
environment, suggesting that baryon fraction is the main driver for bar instabilities. 



expect the bar fraction to be higher in early type spirals, 
that have a higher baryon fraction. This implies that a 
significant fraction of the bulges of early type galaxies 
are pseudo-bulges. 



The morphology-density relation (Dressier 1980) can 
be explained by the notion that the cluster environment 
is cre ating SO's from the infallin g disc population. In- 
deed, iGraham fc Worleyl ( |2008[ ) find that the bulge to 
disc ratios of SU's is similar to that of early type galax- 
ies. One might therefore expect the bar fraction to be the 
same in SOs and early type spirals, however averaged over 
the entire population it is significantly lower (25% versus 
50% respectively). We note that the bar fraction in SO's 
and early type discs with stellar masses above is 



similar 50%), but this drops to less than 10% in the 
least massive SO's. This supports a harassment scenario 
for the formation of the SO population. Gravitational 
encounters between galaxies and with the global cluster 
potential thicken the disks of massive early type spirals 
by an amount that is sufficient to suppress spiral patterns 
(Moore et al 1999). For lower mass disks, the heating is 
more effective and will eventually erase the signatures 
of a preexisting bar. Numerical simulations also indicate 
that infalling late type disks will undergo an environmen- 
tally driven bar instability, however this phase is short 
lived with the bar experiencing subsequent heating until 
it becomes a dE/dSph. 
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